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Abstract: Spectral distortion of the cosmic microwave background provides a unique op- 
portunity to probe primordial perturbations on very small scales by performing large-scale 
measurements. We discuss in a systematic and pedagogic way all the relevant physical phe- 
nomena involved in the production and evolution of the /i-type spectral distortion. Our main 
results agree with previous estimates (in particular we show that a recently found factor of 
3/4 arises from relativistic corrections to the wave energy). We also discuss several sub lead- 
ing corrections such as adiabatic cooling and the effects of bulk viscosity, baryon loading and 
photon heat conduction. Finally we calculate the transfer function for //-distortions between 
the end of the //-era and now. 
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Table 1: Symbols used in the paper. 



1. Introduction 

The cosmic microwave background (CMB) and large-scale structures have allowed us to probe 
primordial perturbations with very high accuracy in the interval of scales 10 -4 < /cMpc < 1. 
Very little is known so far for scales that are parametrically smaller k 3> Mpc -1 . The reason 
is that this information has been a erased by diffusion damping in the CMB temperature 
anisotropy and is swamped by nonlinear gravitational effects in large-scale structures. For- 
tunately the frequency dependence of the CMB spectrum does carry primordial information 
about very small scales. Around the time in which the COBE data became available, a consid- 
erable effort (see e.g. [||, ||, [5], ||]) was devote to the study of the evolution of the CMB spectrum 
and to the search for potential sources of distortion. This effort has decreased considerably 
since the COBE collaboration published data compatible with a blackbody spectrum [pi, im- 
posing strong bounds on the distortion parameters, \y\ < 1.5 x 10~ 5 and \fj,\ < 9 X 10~ 5 . No 
other experiment has since been able to independently improve on COBE/FIRAS constraints. 
Nevertheless the experiment TRIS ||, combining its data with the one of COBE/FIRAS was 
able to slightly tighten the bound on ji, reporting |/x| < 6 x 1CP 5 . Another notable experimen- 
tal effort was ARCADE Q, which detected an anomalous raise in the low- frequency spectrum. 
The subject has enjoyed a recent revival 20, 13, |ll| also thanks to the prospects of an 



observational detection with future experiments such as PIXIE Traditionally, the com- 
putations has been performed at the level of the kinetic theory, using Boltzmann transport 
equation. Given the importance of spectral distortion as probe of primordial perturbations, it 
is useful to have different but equivalent physical descriptions of the relevant processes. The 
fluid limit of Boltzmann transport equation is a very intuitive one and has proven to be a good 
aid to physical intuition in the study of the CMB (for some analytical work emphasizing the 
fluid dynamic approach see |l6|] and references therein). The goal of this paper is to provide 



- 3 - 



a pedagogical, thorough and systematic description of the creation and evolution of //-type 
spectral distortion using fluid dynamics. 

A picture of the phenomena that we intend to describe is the following. At early times, 
z <C Zi = 2 x 10 6 , thermodynamic equilibrium guarantees that photons are very well described 
by a blackbody spectrum. Between Zi and zj = 5 x 10 4 , photon number changing processes are 
frozen out due to the cosmic expansion, but equilibrium is still ensured by elastic Compton 
scattering. During this era, the dissipation of acoustic modes due to viscosity leads to an 
increase of the entropy and the energy of the photon-baryon-electron plasma. Since this 
increase is not accompanied by an appropriate increase in the number of photons, both the 
average entropy and energy per photon grow. In thermodynamic equilibrium, these are given 
by 

tL = T^ [l + p(A n -A p ) + 0^ 2 )] 
ra 7 b R 

~ 2.7 T [l + 0.26/2 + (p 2 )] , 

^ = ~[l + v(A n -A s ) + 0(^)] (1.1) 

~ 3.6 [1 + O.53/t + (/t 2 )] , 

where or, 6r, A p , A n , A s are all numerical constants of order one define in table [l]. Therefore 
an increase in the entropy per photon is tantamount to the creation of a photon chemical 
potential 1 /i. In other words, \x is a direct measurement of how much the rest-frame entropy 
per photon deviates from the fixed blackbody value. An analogous statement is valid also 
for the energy per photon but provided one specifies that T is held fixed or some equivalent 
condition. It is for this reason that s 7 /n 7 is a more useful indicator of /i-distortion than 
Pj/n 7 . 

After Zf, thermodynamic equilibrium fails to be reached because the plasma temperature 
has dropped well below the electron mass. From this moment onward, the creation of /i-type 
distortion is halted (while y-distortion can still be generated but the distortion present 
at the end of the p-eia, zt, does survive and evolve in a non-trivial way all the way to us. On 
large scales, this evolution consists simply of the projection of the inhomogeneities present 
at Zf onto the multipoles of the late time CMB sky. On smaller scales, I > 1000, diffusion 
damping and the thickness of the last scattering surface erase any inhomogeneity. This last 
result is relevant even if we were not to observe these small scales. To understand the reason 
let us assume a Gaussian universe with just adiabatic initial conditions. Then fluctuations of 
the /i generated by the dissipation of acoustic modes ((/t) ~ few 10 -8 ) on large scales are just 
white noise. Hence the /i-power spectrum on large scales, C^ 1000 , is linearly proportional to 
the largest dissipation scale because of an average along the line of sight. This was derived 
in [11 1, where the largest dissipation scale was called k s , and will be discussed in section |j. 
The evolution after the /i-era erases inhomogeneities in /t also on scales that are larger than 



x We use the symbol /i for the dimensionless size of spectral distortion of the fj, type. This is related to the 
standard (dimensionful) chemical potential /t t h used in thermodynamics by fi = —fith/T. 
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the dissipation scale at the end of the //-era due to the finite thickness of the last scattering 
surface and further diffusion processes. This suppresses the large scale-power spectrum of 
\i even more. Hence, in a Gaussian and "adiabatic" universe, there is practically no hope 
to extract any information from the \x multipoles. On the other hand, in the presence of 
primordial non-Gaussianity 2 things are different fll|| . [i is modulated by large-scale potential 
fluctuations which also produce the large scale temperature anisotropics. Thus its correlation 
with 5T gives a direct measurement of primordial non-Gaussianity. Prom this point of view, 
the Gaussian power spectrum discussed above is just some source of noise, typically referred 
to as cosmic variance. The fact that this cosmic variance is negligibly tiny for [x anisotropics 
(generated by the dissipation of acoustic modes) can be of great advantage in constraining or 
detecting deviations from Gaussianity at very small scales. 



1.1 Summary 

Let us give a quick summary of our main results and of the organization of the paper. In 
section ^ we set the stage with a pedagogical review of the hydrodynamic description of 
the universe much before the //-era. All the results derived there are already present in the 
literature. In section || we adopt the same tools to provide a detailed descriptions of the //-era 
and the evolution of the /i-type spectral distortion. A lightning derivation of the main result 
goes as follows 3 . We expand in perturbations around an FLRW universe containing just a 
photon-electron-baryon plasma 4 . The evolution equations are the conservation of the total 
energy momentum tensor T-ji = 0, the conservation of baryon number Njf = and, during 
the /i-era, the conservation of photon number Ny ilM = 0. At linear order in perturbations these 
equations describe the propagation of waves. Including dissipative corrections, the waves get 
damped once they reach some small scale qo- Perturbations to the photon density <5 7 and 
velocity potential u (so that the irrotational fluid velocity is given by Ui = d{u) are of the 
form 



1 + w 



qu(t,q) ^S 7 (t,q) ~ A cos (J ^c s dt'^j e~ $ Tdt ' , (1.2) 



where A is some amplitude determined by the superhorizon primordial perturbations, c 2 . ~ 
w ~ 1/3 are the speed of sound of the fluid and the background expansion index, respectively, 
and r is the dissipation rate given explicitly in ( |2.69f) . There are corrections of order the 



2 Also in the presence of isocurvature modes things can be more interesting, but we do not discuss this case 
in the following 

3 This derivation uses rest-frame quantities such as n 7 and p 7 . A parallel derivation using cosmological- 
frame quantities such as N® and T 00 is given in appendix |^. The two results are identical for /i, since this 
quantity does not change as we move to a different reference frame. The results for temperature fluctuations, 
on the other hand, are different but related by the transformation in ( F.16| ). 



4 Dark matter is indeed a small corrections since the /t-era happens entirely during radiation domination. 
The presence of neutrinos, on the other hand, leads to a 10% error that we fix in the more rigorous treatment 
in the rest of the paper. 



- 5 - 



baryon loading R, which we neglect in this summary, but will be accounted for in the following 
sections. In any case, these corrections are small, since during the /i-era R -C 10~ 2 . 



To study the evolution of [i we first notice that according to (LI) 

4a R 



t 



n- 



1.9/2. 



(1.3) 



Since according to COBE/FIRAS Q /x <^ 1, here and in the rest of the paper we work at 
linear order in /i. To compute s 7 and n 7 we use the conservation of photon number and the 
evolution of the entropy density 



ATM 

7;m 



I SW 



0, 

-AT' 1 "!* 



AN v d uf x , 



(1.4) 
(1.5) 



where the latter equation is derived in appendix |^. Here is the fluid velocity of energy 
transport (see e.g. section |2|) while AA^ 4 and As^ are viscous corrections that account for the 
fact that, for an imperfect fluid, the velocities of the transport of photon number and entropy 
can be different from u^. As we show in appendix O, up to terms of order R, ANtf = 0, so 
we can focus on AT^ V . Its general form at leading order in derivatives is given in (2.51) and 
it is proportional to the shear viscosity r), which for the photon-electron-baryon fluid takes 
the value [25| r\ = (16/45)t 7/ o 7 , where t~ is the photon mean free path. Solving 5 (|1.4 



and 



1.5) for s'ry and ri 7 one finds 



7U 



nT M ' 



(1.6) 



Since the left hand side is proportional to /I, this implies that /i-distortion is created only 
when the viscous corrections are important (again this is true only up to corrections of order 
n b/ n -y)- As discusses in subsection 2.2.1, AT^ V starts at linear order in perturbations and 
0. Then the leading term in the right hand side of (|1.6|) is second order 



ArpOfM 



in the small amplitude of perturbations. The explicit computation using ( |3.14| ), ( 2.61 ) and 
(P| gives 



15(^4n A s 



(d lUj ) (sv + av) - ^ (V 2 



(1.7) 



If we take the ensemble average, focusing on the average /i-distortion, and use the relation 



between u and 5 7 ([[j) and the explicit expression for T ( 2.6E ), this expression agrees with 
the result of fl2l 



(A) = -lAd t ((E w /p)) 
l + w ' 



(1.8) 
(1.9) 



3 This is most easily done in the rest frame of the fluid where <9 M u° = u l = and u° — 1. 
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where E w is the energy of a wave, () the expectation value over the ensemble of all histories 
and ()p is the average over a period. Notice however that it is important to include the correct 
relativistic factor (1 + w)^ 1 ~ 3/4 (derived in appendix |E|) in the formula for E w . Therefore 
the derivation of the correct numerical factor (first noticed in j^0|) can be obtained entirely 
within the realm of fluid dynamics, without any reference to the microscopic theory. 

In the more detailed description of section |3| we identify other subleading contributions 
to /i-distortion: 



• The baryon temperature leads to adiabatic cooling (subsections |2jj and 3J) [24] and a 
non- vanishing but small bulk viscosity (subsection |2.2j ). 

• The baryon energy density (mostly rest mass energy) leads to a small damping of the 
amplitude of acoustic waves (subsections |2.1,2 , 2.2.3| and 3.3.3| ) and is small during the 
//-era, but could be relevant for the y-era. 

• Similarly to what happens for baryons, the conservation of photon number is corrected 
by heat conduction (derived in appendix [D]) . 

After the end of the /i-eia, the dissipation of acoustic modes cease to produce spectral 
distortion of the /U-type. Nevertheless, in order to make predictions for late-time observables, 
we need to evolve the signal up to the last scattering surface and then account for the free 
streaming of photons all the way to us. We do this in section ||. The final result ( 4.14| ), is 



simple to understand, /x-distortion, defined through the dimensionless chemical potential /j, 
is not affected by free streaming or gravitational potentials. Hence for I < 1000 the transfer 
function is just the simple projection of the /x monopole at the last scattering surface onto 
various multipoles as in ( 4.17| ). For / > 1000 two effects become important: the finite thickness 



of the last surface ( 4.14j ) and diffusion damping (4.5). The outcome is an exponential damping 



of /i anisotropics for I > 1000, and a suppression of its power spectrum, as discusses in section 



2. Before the /i-era 

To set the stage, in this section we describe the evolution of the universe at early times, 
z 3> Zi = 2 x 10 6 , when double Compton and Bremssstrahlung interactions of photons with 
the electron-baryon plasma were very efficient. At that time almost perfect thermodynamical 
equilibrium was maintained and the photon occupation number was very close to a black 
body spectrum. Since photons can be efficiently produced, the only conserved number in 
this black-body era is the baryon number. After introducing the relevant hydro dynamical 
equations, we derive the evolution of primordial adiabatic perturbations at linear order and 
present those second order result that are relevant for the production of spectral distortion. 

Because of the frequent interactions, photons electrons and baryons combine into a single 
fluid, whose energy momentum tensor is covariantly conserved 

T.^ = 0. (2.1) 
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These equation have to be supplemented by the conservation equations of any conserved 
charge in the system. After baryogenesis, when the temperature has dropped well below a 
GeV, the baryon number is conserved. Later, after electron positron annihilation, z ~ 10 9 , 
also the electron number becomes effectively conserved. However ti a very good approximation 
matter behaves as a single fluid, so it suffices to follow baryon number conservation 

= > ( 2 - 2 ) 

where is the baryon- number current, defined such that in the rest frame of the fluid N° 
corresponds to the baryon-number density. Finally we need to consider Einstein Equations 

GV = , (2.3) 

where M v = (8ttGn) is the reduced Planck mass. 

2.1 Inviscid fluid 

The energy-momentum tensor for a perfect fluid is diagonal and isotropic in the rest frame 
of the fluid, hence in any frame it is given by 

= (p + p) u^Uy + g^p , (2.4) 

where g^ v is the metric, is the four- velocity of a fluid element with u^u^ = — 1 , the energy 
density is the one measured in the rest frame of the fluid p = u^UyT^ and for a barotropic 
fluid, the equation of state gives the pressure p as function of p. As mentioned previously, 
in the presence of a conserved quantity, such as baryon number, there is a Lorentz frame in 
which = {n,0}, with n the number density in the rest frame. For a perfect fluid, the 
energy-momentum tensor is diagonal in this same Lorentz frame, hence in any frame we can 
write 

iV" = nu^ , (2.5) 

so that n = —u^N^ 1 . For baryons the pressure, number density and energy density are related 
by 6 

3 

p b = n b T , p b = -mn b + Tn b . (2.6) 

Here m ~ GeV is the average mass per baryon and T is the baryons temperature, which until 
late times, z ~ 200, differs little from the photon temperature. We will neglect this difference 
and use T for the whole photon-baryon-electron fluid. For photons, the pressure is given by 

p 7 = p 7 /3, (2.7) 



6 We set the Boltzmann constant to one, fcs = 1, i.e. measure temperature in units of energy. 
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independently on whether the photon number is conserved or not. For a black body spectrum 



Pl = a R T A , , (2.8) 
with or = 7r 2 /15. For the photon-baryon-electron fluid we then take 

p = p 7 + mn b + Tn h , p = \-p^ + n b T . (2.9) 

Notice that baryon temperature is typically a small correction. On the one hand, as long as 
we are interested in the time much after big bang nucleosynthesis (BBN), T <C -1 MeV, and 



so we can neglect the temperature dependent term in the baryon energy density (2.6), up to 
corrections of the order T/m <C 10 -4 . On the other hand 

r(t) ^PJL^ = ^^gl ^ .37 x^,2.2x 1(T 10 , (2.10) 

p-y p 7 n 7 7r 4 n 7 

with C the Riemann zeta function and n 7 the number density of photons. Despite being time 
dependent, the number of photons per baryon, which we parameterize by r, has changed very 
little since electron-positron annihilation. BBN and CMB bounds agree, within the error 
bars, on the value ri&/n 7 ~ 6 x 10~ 10 . 

Because its effects are suppressed by these small factors, baryon temperature can play an 
important role only when all other contributions cancel precisely. This does indeed happen for 
the spectral distortion generated by the homogeneous adiabatic evolution during the //-era, 
as we will review in sections |2.1.1| and |3.3| . In all other sections but these two, we can safely 
neglect terms proportional to the baryon temperature. 

We now proceed to solve ([O]) and fl2.2|) in a perturbative expansion around a flat FLRW 
universe. It will be convenient to introduce the notation 7 



Pb^ — Pb,-y 



1 + <^,7 + C 



(2.12) 



Ui = Ui + Ui + uf> , ( 2 -13) 



"6,7 = n W 



1 + 5n bn + dnf^ 



(2.14) 



so that any 5 and 5n are dimensionless. A bar indicates zeroth order homogeneous (but 
time dependent) quantities and to simplify the notation we omit label "(1)" in the first order 
perturbations. Finally notice that due to the tight coupling between photons and electron- 
baryons, there is a single fluid velocity U{. Since vector perturbations decay in and FLRW, 
we neglect vorticity and introduce the velocity potential u by 

U{ = d{U . (2-15) 



7 Some authors (e.g. E|) define 5 7 = 5p 7 /(p 7 + p 7 ). Our choice of 5 is related to that choice by 



= (1 + R)4 S 

7 3 + 4R 7 v ' 



where R = 3pb/(4p 7 ) is the baryon loading. 
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2.1.1 Zeroth order 



We consider a homogeneous and isotropic background, so that p = p(t) and Uj = 0, and the 
unperturbed metric is the one of an FLRW universe 



ds z 



-dt 2 + a 2 dx 2 , 



with a(t) the scale factor. The Christoffel symbols are found to be 



1 jo 



0j 



Hsu , r: 



a 2 H5, 



(2.16) 



(2.17) 



where H(t) = a /a is the time-dependent Hubble parameter. Einstein equations then reduce 
to the Friedamann equation 



3M 2 H 2 = p~/ + p b + p c + Pu , 



(2.18) 



where p c>v are the dark matter and neutrino contributions, respectively. Using ( |2.17| ), the 
number (^^) and energy conservations ( |2.lD lead to 



n b 
Pi 



AHp y 



1 + -r + O (r 2 ) 



(2.19) 
(2.20) 



where r was defined in (2.10). Since there are about 10 9 photons per baryon, the baryon 
pressure can be neglected up to corrections of order r ~ 2.2 x 10 -10 . In (2.20) we see that, 



the effect of the baryon temperature is to accelerate the dilution of the radiation energy 
density. This can be understood (see e.g. [^qj ) by recalling that if the photons and the baryon 
were not interacting, their temperature would redshift as a -1 and a~ 2 , respectively. Once 
the interaction is taken into account, the common temperature scales with a power of a in 
between —1 and —2. Since the are so many more photons than baryons, the scaling in actually 
very close to —1, as can be seen by rewriting ( 2. 20] ) using (2.8) 



-HT 



l + -r + (r 2 ) 



(2.21) 



2.1.2 First order 

For the study of perturbations we choose to work in synchronous gauge 8 

„2 i.2 , 2.J. 



and define 



ds z = -dt z + a z dx l dx 3 [(1 + A) 5ij + B,y] , 



(2.22) 



(2.23) 



3 We use the notation of Jl(|, except for the definition of the density perturbation as explained in footnote 
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With this choice, at linear order the time component of the velocity vanishes 5u° = 5goo/2 = 
and no metric perturbation appears in Uj = a 2 u % . Also the only non- vanishing metric 
perturbations are 5g lJ = —a~ 4 5gij, so indices of metric perturbations are not raised and 
lowered by the background metric. 

Let us now expand equations (|2.1|) and (|2.2|) at linear order in perturbations 



d t T$ + 3HT§ + a 2 HT i { [ ) 






1 (0) 01 pLCt 1 


(2.24) 


d t T^ + 5ifTg 


+ ^ 




± (o) oi M a ' 


(2.25) 




+ <K ) 


= 0. 




(2.26) 



Neglecting the baryon temperature, i.e. terms of order r, the energy-momentum tensor at 
linear order is given by 



rpOi 

1 W 

rrAj 
i (l) 



P 7 ( <5 7 + -RS b 



-/9 7 a 



-2 



(<%5 7 



-2 



(2.27) 

(2.28) 
(2.29) 



and 5T can be found in textbooks, e.g. [16]. Going to Fourier space, ( [2.24 -2.26) become 



: (R + l)u = -(R+lty, 



(R+ l)u - Hu+ -Sj 
on;, 



(2.30) 
(2.31) 
(2.32) 



where R = 3pf ) /(4p 1 ) is the baryon loading, which to leading order is just proportional to the 
Cinstein e 



scale factor a. While Einstein equations give [16] 

1 



(2p 7 (5 7 + p b 5 b + p c 5 c + 2p u 8 v ) . 



(2.33) 



One can solve (|2.32 ) for 5n b = 5 b , (|2.30| ) for u, ( [2.31 ) for u and plug them into the time 
derivative of (|2.30| ) . The result is 



■i rr 1 + 2R 

7 l + R 



(X O 



(2.34) 



where we have introduced the speed of sound in the photon-baryon-electron fluid c 2 = 
[3 (l + R)]' 1 . Analytical solutions of (^3§ and fl2~3!p can be quite involved, but as long 
as we are interested in the evolution of perturbations well inside the horizon things are much 
simpler. During radiation domination, the gravitational perturbation tp decays inside of the 
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horizon as q 2 /(aH) 2 due to the large pressure of photons, so we can drop all ip terms in fl2.34|) . 



Then (2.34) can be easily solved by the WKB approximation. 



We still need to find the correct normalization of <5 7 fluctuations inside the horizon in 
terms of primordial superhorizon fluctuations. Here we consider only adiabatic initial con- 



ditions, specified in terms of the gauge invariant variable 7Z defined in (B^), which in the 
simplest cosmological scenarios is conserved outside the horizon. If we neglect dark matter 
and neutrino, then the universe is filled with a single fluid and can be described using a 
P{X) Lagrangian as in []. The derivation of the initial condition, whose details are collected 
in appendix [B], then proceeds as follows. One works with the gauge invariant variable 1Z 
whose action ( |B.1| ) is well know. Its equations of motion can be easily solved both inside and 
outside the horizon, during radiation domination. Then well inside the horizon one goes to 
synchronous gauge in which 7Z and the velocity potential u are directly related. Finally one 
can relate this solution for u to a solution for <5 7 using ( |2.30| ). The final result is 

8 -> {t) = w ™ (/' W) c ' {t ' )dt ') + (I) ' (2 ' 35) 

with 1Z° the asymptotic value of 1Z at early times, outside of the horizon, which is determined 
by the microphysics of the early universe, e.g. by inflation. The factor (1 + it!) -1 / 4 is best 



understood 9 by computing the leading order WKB solution of ( 2.34 ). 

The above derivation does not capture the dark matter and neutrino contributions. While 
the first component is indeed very small during radiation domination, the second leads to a 



10% correction. The numerical analysis of [18] gives 

= ~ /4 (1 - 0.268^) cos (^J -S-c s (t')dA . (2.36) 

Here and in ( |2.35| ) we have re-absorbed a neutrino dependent phase into the lower bound of 
the integral. The quantities we are interested in contain the average of oscillations over a 
period, hence the phase will be irrelevant for us. Here we have defined ratio of neutrino to 
photon energy density 

R -Tr N »\{rT- (2 ' 37) 

For the standard model with three neutrinos, N v = 3.04 and one finds R u = 0.41. Combining 
( p0| ) and ( p2|) , one finds 



\u = \ <j 7 + V • (2.38) 



9 One might be concerned that terms of order R were dropped in computing the normalization of <5 7 in 
terms of the conserved quantity 1Z. Notice however that R grows with time, so it makes sense to neglect terms 
of order i? C ari y and keep terms of order -Riato 2> -Reariy 
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If we focus on modes well inside the horizon and neglect tp, to leading order in aH/q we get 



d 3 q \TZ°(q)\ : 



3c s (1 - 0.268^) sin 



(2vr) 3 (1 + R) 1 / 2 
where by statistical isotropy and homogeneity 

{K (q)K°tfr) = (2vr) 3 <5 3 (?+ , 
For future reference we also compute 

d 3 q \K°{q)\ 2 



-Co dt' 



\K\q)\\ 



(S 



(2vr) 3 (1 + fl)i/2 



4(1 - 0.268i2„) cos 



-c, dt' 



(2.39) 



(2.40) 



(2.41) 



A few comments are in order. First, taken at face value, given a scale invariant spectrum of 
primordial perturbations both integrals in ( 2,39j ) and ( p. 41 ) are logarithmically divergent (in 
the UV and IR). We will see that for the physically relevant quantity, i.e. fi, the divergences 
are naturally cut off. Second, notice that the velocity potential and density perturbations 
oscillate with a phase difference of tt/2. Hence, if we denote by () p the time average over one 
oscillation (while () is still the quantum expectation value) we find 



{{u l Ui)) p 



(2.42) 



Finally, the amplitude of density and velocity waves, i.e. ((u l Ui)) p and {{S 2 )) p , depends very 
weakly on time. Deviations from a constant amplitude come at order R, which is very small 



during the /U-era. We will come back to this point at the end of subsection 2.2.3 



2.1.3 Second order 

Let us move on to second order perturbation theory. A thorough study of second-order 
perturbation theory is beyond the scope of this work. We will instead concentrate on the 
results relevant for the production of spectral distortion, i.e. we want to study how the 
dynamics of short waves feeds back at quadratic order into the homogeneous evolution. First 
of all we will neglect gravity perturbations. The reason is that we will be interested in the 
physics of perturbations that are well inside the horizon where gravitational potentials are 
negligible. A great simplification is achieved by studying the evolution of the expectation 
values () of the various perturbations, rather than the perturbations themselves. Assuming 
statistical homogeneity and isotropy one can then discard all derivative terms. 

Let us take the expectation value 10 of the conservation of energy {T?^.^) = and particle 
number (Nff) = 0. Spatial derivatives can be taken out of the expectation value. Then 
invoking statistical isotropy they must all vanish. One is left with: 



((3H + d t )T^+a 2 HT^) 
((3H + d t )N° 2) ) 



0. 
0. 



(2.43) 
(2.44) 



'This can be thought of as a quantum/statistical expectation value or as an average over space. 
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where 



1 (2) 



J (2) 



P- 



r^ 7 (2) + (1 + R)uiU l + RS b (2) 



2 <%7* 7 (2) + + 



iV ( Q 2) = n b 5n b(2) + -u l 



(2.45) 
(2.46) 
(2.47) 



and we used it?™ = u i Ui/2. We focus on the ensemble average and use statistical homogeneity 



4 (2) 

drop all total space derivatives. Then 

,3 , 



(7*7(2) + (1 + ("it**) + i2«5 6 (2) + RHu lU l 



dt ( ($.(2) + ^i^> 



0. 
0. 



(2.48) 



We differ a discussion of this result to subsection 2.2.3 , where we will be able to include the 
effects of diffusion damping. 

2.2 Viscous fluid 

The perfect fluid provides just the zeroth order behavior in an expansion in small momenta. 
Every realistic fluid deviates from this idealized case and deviations can be systematically 
parameterized in terms of viscous coefficients. For our purposes it will suffice to consider 
just the first order corrections in the derivative expansion. In general one will have (see 
e.g. §§0) 

where T^ v and are still given by the perfect fluid formulae fl2.4| ) and ( |2 . 5|) . One can 
still define p and n as the energy and number density in the rest frame defined by u^, 
i.e. p = u^u u (T^ u + AT^ V ) and n = -u^ (iV M + AN*). This implies that in any frame 



UaU u AT^ = , UuAJV = . 



(2.50) 



Then one can still redefine such that either T°* + AT * or N l + AN 1 vanish in the rest 
frame defined by u^. We choose the first option so that in any frame u (J AT ' 1 = 0. In words, 
this means that is the velocity of energy transport. Because of this choice the coarse 
grained fluid velocity of particles whose number is conserved, e.g. baryons, is not the velocity 
u^. Hence their conservation has to be corrected by a non- vanishing AN 1 - 1 . The generic 
corrections can be parameterized by three coefficients according to 



AN^ 



-X 



nT 
P + P 



(2.51) 
(2.52) 



- 14 - 



where rj, £ and x are the shear viscosity, bulk viscosity and heat conduction respectively. 



The second law of thermodynamics requires them to be all non-negative. In (2.52) we have 
defined the dimensionless chemical potential p by dividing the standard chemical potential 
by T. This definition agrees with the common use in the literature on spectral distortion and 
is quite convenient for photons, as we will see. 

The viscous coefficients rj, ( and Xb, where Xb is the heat conductivity appearing in the 
conservation of baryon number have been computed 11 in [fTi 



16 

45 



P'jt^ , XbT r, Ay^7 J C 4/9-yt,y 



dp 

d pJn. 



(2.53) 



Here i 7 = l/(<7Tn e ) is the photon mean free path, with n e the density of free electrons of mass 
m e and charge e and ot = 87re 4 / (3m , 2 ) ~ 66 fm 2 is the Thomson cross section. In appendix 
|d] we rederive Xb as well as the analogous heat conduction for photons Xj which is relevant 
during the ^-era, when the number of (not-so-soft) photons is conserved. The bulk viscosity 
is extremely small but non-vanishing because of the baryon temperature. From ( |2.53| ) one 
finds 



( = p 7 t 7 - + 0(r 3 ) ^T Xb , V- 



(2.54) 



2.2.1 Zeroth Order 



As in the perfect fluid case, we choose a homogeneous and isotropic background with p = p 



and Ui 



0. By the definition of p and one must have AT^ = 



By isotropy 



and the only non vanishing component is then the spatial diagonal part 



AT ij = -£ a -\.sH . 



(2.55) 



We see that even at the zeroth order, entropy could be created 17] in the presence of a sizable 
bulk viscosity. Since Q for the photon-baryon-electron plasma is so small (see ( p. 54 )), we can 
safely neglect this effect, up corrections of order r 2 ~ 10 -19 . In conclusion the background 
evolution is unaltered from our perfect fluid discussion of subsection |2.1.1| . 

Before proceeding, let us see which other constraints can be derived on the viscous con- 
tribution. By using 



u 







1 . 



(o) — x ' u (l) — u ' "(0) 
= u^AN^ at subsequent orders in perturbation theory, one finds 



o 



0. 



0. 



0. 



0. 



and Ufl AT^ 



(2.56) 



AT ( °i) 



AiV { ° 2) 



0. 



(2.57) 



L The shear viscosity in jl7| should be corrected for the contribution of photons polarization JlC 
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2.2.2 First order 



At linear order in perturbations, we want to compute AT^ U . These viscous terms in general 
become relevant only at some dissipation scale which is approximately given by the geometric 
mean between the Hubble and mean free path scales. Physically this corresponds to dis- 
tance that a typical microscopic particle covers due to Brownian motion in a time H' 1 f25fl. 



Given that t 7 <C H at any time before decoupling, the dissipation scale is well inside the 
horizon. During radiation domination the gravitational potentials decay inside the horizon 
as q 2 /(aH) 2 , so we safely neglect them in computing AT and AN. In order to use ( |2.52|) we 
need to compute the chemical potential for baryons. This is straightforward once we have 
any one thermodynamical potential in terms of its natural variables. A simple choice is to 
compute the entropy from the microcanonical potential (Sackur- Tetrode equation). Ignoring 
very small quantum effects one finds 12 



M6,th 



dS(U, V, N) 
dN 



uy 



m 



11 



(mT)3/ 2 



+ -log(2vr) 



(2.58) 



Given that after BBN T <C m, one can neglect the logarithmic contributions. Using 2.57 and 
neglecting baryon temperature, i.e. corrections of order r, we find 



ATP U U = a'HATfis = , 



(i) 



AT™ = djAT^ 



0; 



-r/a ( djU 1 + diU J — -Sijdiu ) — (a 5ijdiu 



AN^ = a^df, (a 3 AN») 

= - * bf R 
~ n " fa (1 + Rf 



(2.59) 
(2.60) 

(2.61) 
(2.62) 

(2.63) 



Now we set £ = 0, use the zeroth order solution T ~ —HT, trade temperature for density 
perturbations AST /T = <5 7 and focus on scalar perturbations u % = a~ 2 Ui = diU. The result is 



AT. 



(2.64) 
(2.65) 



Because of ( [2.59| ) the energy conservation equation ( 2.24 ) is unaltered, while the momentum 



12 Notice that the relativistic chemical potential is related to the non-relativistic one by fith = Mnon-rci + m. 
This is because in the non-relativistic description the energy does not contain the rest mass of the particle. 



- 16 - 



and number conservation in Fourier space become 

2 

^(R + l)u = -(R+l)iP, 
1 

4' 



3 ■ 

-(L + R5b 

4 a' 



(R + l)u - Hu + -5 y 



-2T v u , 



R ( 5ni, — 



-Rip + 6r 



X 



Hu 



where, with some hindsight, we have introduce the well-known damping rates 



r - 1 — 



R? 



a 2 6(1 + i?) s 



7 a 2 (l + R) 45 ' 



(2.66) 
(2.67) 
(2.68) 

(2.69) 



As in the perfect fluid case, we can solve for 5rib = 5b, u and u and find a second order 
ordinary differential equation for <L 



5 7 + J 7 



h . „ + 2 (r x + r,,) + o I 



o 2 

+5 — c 2 



2 " X 

fHT x a 2 



+ 



0. 



(2.70) 



The two terms that we have not written down explicitly are small and can be neglected. The 
first is suppressed by (aH/q) 2 which is small inside the horizon, the second is suppressed by 
HT x a 2 /q 2 ~ Ht^ which is small at any time before Hydrogen recombination. The solution of 
this equation can be found with the WKB approximation. At leading order, the real part of 
the time dependent frequency is unchanged with respect to the perfect fluid case, while the 
imaginary part has a new contribution proportional to the sum of the dissipation rates. The 
solution normalized as discussed in section [2.1.2| is 



SJt) 



( 



ATI ( f* 
— (1 - 0.268^) cos / 



(2.71) 



Another convenient way to re-write the damping is by introducing the damping scale 



Id 



dz(l + z) 



R 2 16 
1 + R + 15 



(2.72) 



6(1 + R)H 7 

so that 13 <5 7 ,u cx exp (—q 2 /q D )- As we did in (|2.38|) , we can use this result to compute u l Ui. 
If we are interested in modes well inside the horizon we can neglect tp and terms suppressed 
by aH/q. As long as the wavelength of the modes is also much longer than the mean free 



path t 7 we can also neglect the heat conduction term in ( [2.68| ) . Then we find 



{u l Ui) 



d 3 q \TZ°\ 2 
(2vr) 3 (1 + R) 1 / 2 

/ r 2 



3c s (1 - 0.268^) sin (J* ^c s dA e "f ( r x+W 



i 2 H 2 ' a 1 



(2.73) 



13 Notice that with this convention |2q ], the damping scale in the power spectrum will have an additional 
factor of \/2- 
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All the discussion at the end of subsection 2.1.2 carries through, in particular fl2.42j) is valid 
also in the presence of dissipation. 



2.2.3 Second order 

In principle viscous corrections can affect the dynamics of second-order perturbations which 
we studied in subsection 2.1.3] in the perfect fluid case. In practice though, one expects viscous 



corrections to become relevant only at scales that are parametrically smaller than the horizon, 
since these are all higher-order in an expansion in derivatives. As we discussed in section |l], 
all of the modes that we can realistically hope to detect, i.e. I < 1400, where well outside 
of the horizon for the whole duration of the //-era. For these modes we can then neglect 
the effect of viscosity, except for taking them into account in the solution of the first-order 
perturbations that source the second-order evolution. Hence, the solution is again the one 
found in subsection 2.1.3 , i.e. ( [2.48 ), where now for uiu 1 we should use ( |2,73| ), Focussing on 



<5 7 ( 2 ) one finds 

(<* 7 (2)) = ~|(2 + RWunj) - ^HR^iU 1 ) . (2.74) 

This equation tells us how the evolution of the background average energy density changes 
due to the dynamics of acoustic waves. There are three types of terms which source (# 7 (2)) and 
hence represent the backreaction of waves on the homogenous evolution: dissipation terms 
and oscillation terms that survive in the limit R — > and terms of order R. 

• Dissipation terms. These are the most important and are present also for R = 0. 
They arise when the time derivative in ( |2.74j ) hits the damping factor in ( 2.73| ), The 



physical picture is that dissipation erases all waves on small scales and homogeneously 
redistribute their energy, hence sourcing (5 7 ( 2 ))- 



Oscillation terms. They arise when the time derivative in ( 2.74 ) hits the sines in fl2.73|) . 



These terms are present also for R = 0. The physical implication is that even if the 
amplitude of the waves were constant, there would be oscillations in (<5 7 (2))- Notice 
that the average energy density (T&y is constant and oscillations arise just because p 
is defined as the energy density in the rest frame rather than the cosmological one. 

• Terms of order R are always small corrections during the /i-era since R increases mono- 
tonically with time reaching its maximum at the end, when R(zf) ~ 0.01. In addition, 



as we will see in subsection \$.2[ all these terms should cancel in the source term for (fi) . 



We will provide more details on the dissipation terms in 2.2.3, when we will discuss the 
generation of spectral distortion. 

3. During the //-era 

As the universe expands, double Compton and Bremssstrahlung interactions become less and 
less efficient. These interactions are the leading processes that change the number of photons, 
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since elastic Compton scattering does not create new photons. After redshift of z% = 2 x 10 6 
double Compton and Bremssstrahlung interactions can efficiently create new photons only at 
very low frequencies. But Compton scattering is not fast enough to redistribute these photons 
to higher frequencies. The end effect is that the number of photons above a certain low- 
frequency threshold becomes effectively conserved, i.e. the //-era begins. We can still describe 
the system using hydrodynamical equations, but the conservation of energy, momentum and 
baryon number needs to be supplemented by another equation for the conservation of the 
number of photons 

(A^ = 0. (3.1) 

During the /i-era Compton scattering is very efficient at exchanging momentum between 
photons and electrons (and eventually baryons). This means that perfect thermodynamic 
equilibrium is reached but with a conserved charge, i.e. the number of photons. By studying 
the Boltzmann equation with just Compton scattering in the collision term (Kompaneets 
equations Q), one finds || |6| that the occupation number of photons is well described by 
the Bose-Einstein distribution with some temperature T and chemical potential fi. Energy 
momentum tensor and number density are given by 

1 f ol 3 p p^p u 1 
T^(x,t) = -3 J (2vr)3 p o e p°/kT(x,t,p)+n(x,t) _ I ' ( 3 - 2 ) 

OM) - ^3 J (2^)3^0 e p°/kT(x,t,p) +f ,(x,t) _ 1 ' ( 3 - 3 ) 

where with this definition /x is dimensionless. We adopt the same sign that is commonly used 
in the literature, i.e. the opposite to the one used in thermodynamics. In the tight coupling 
regime the mean free path is so short that an observer in the rest frame detects an isotropic 
distribution, i.e. T and \x do not depend on p, perturbations have only a monopole. Also T 
by definition does not depend on \p\ and /i can found to be 



with Ppeak ~ !0~ 2 kT, so that it is very close to a constant for middle and high frequencies. 
In the following we will simply use [i to indicate this constant value. Then the rest-frame 
photon number density is 



"^rest 



d 3 p 1 



(2tt) 2 e P°/T+M - 1 
[1 + 3G r + 36^(2) + 39^ - A nf i] . (3.5) 



where A n = -vr 2 /[6C(3)], p°p° = p^aT 2 = p 2 a~ 2 , n 7 = T 3 b R with b R = 2C(3)/vr 2 , and we 
defined the rest-frame temperature perturbation by @ r = 5T/T. In the second line we have 
expanded to linear order in [i and quadratic in 0. Terms quadratic in [i are very small and 
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we neglect them. A similar computation leads to 



rpOO 

rest 



Pi 



.-3 



d 3 p 



P 



(2vr) 2 e P°/^+M _ 1 
p 7 [1 + 46 r + 4e r)(2) + 69^ 



(3.6) 



with = 90C(3)/7r 4 , |0 7 = orT 4 . As the universe keeps expanding the temperature as well as 
the density of electrons decreases. After redshift of about Zf = 5 x 10 4 the interactions become 
so rare and the typical momentum exchanged so small that elastic Compton scattering is not 
efficient enough in maintaining thermodynamical equilibrium 14 . This time signals the end of 
the /i-era. At any time after Zf ~ 5 x 10 4 , any perturbation to the system cannot be efficiently 
thermalized. The type of distortion arising at later times is of the y-type and has a different 
frequency dependence from the distortion generated during the //-era. 



3.1 Inviscid fluid 



In this section we will repeat the computations performed in section 2.1 supplementing the 



equations of motion with an equation for the conservation of the photon number. The main 



result is that fi is sourced only by terms of order nb/n 7 , which was estimated in (2.10). Let 
us to start considering the idealized case in which photons behave as a perfect fluid. We can 
write the photon number current as 

Nt; = , (3.7) 

where n 7 is the photon number density in the rest frame of the fluid given in ([T^). For an 
imperfect fluid, this expression is corrected by a heat conduction term. This will be discussed 
in section [3]3] and appendix [D| and turns out to be a negligible effect. 

We study first the evolution of the zeroth-order homogeneous background. As long as we 
neglect viscous corrections, the entropy is conserved and the system evolves along an adiabat. 
This is true both at the homogeneous and inhomogeneous level, therefore the zeroth order 
results will be sufficient to understand the behavior at higher orders in perturbation theory. 
The number of photons is conserved in a comoving volume 

rt 7 + 3fln 7 = . (3.8) 

Adding this to the equations of subsection |2.1.1| , and using (3J3) and (|3.6|) to rewrite n 7 and 



14 This is sometimes referred to as kinetic equilibrium to stress the fact that interactions isotropy photons 
but do not change the amplitude of their momentum. 
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p 7 in terms of 5T and /i, on finds 



rib = —3Hrib , (3-9) 
f 1 = -Hf 



1 + - r - + O (r 2 ) 

+ 8 tt 6 - 405C(3) 2 ^ { ' 



-HT I 1 + -r2.5 ) , (3.10) 



^ 4 vr6 - 405C(3) 2 + 1 j ' 
~ -2.1rH. (3.11) 

So we find that a homogeneous negative /U-distortion is created by the adiabatic expansion. 



This effect to was first noticed in [24] and it is sometimes referred to as production of fi- 
distortion by adiabatic cooling of electrons. Intuitively, electrons and baryons would tend to 
cool down faster than the photons. Because of the continuous interactions, the electrons and 
baryons extract energy out of the photons. This results in a negative /^-distortion. As the 
name suggests, no entropy is created in this process. Notice also that the amount by which 
the temperature evolution deviates from a~ l is different from what we found before the //-era 



(2.21). The typical amount of /i-distortion produced by this effect is of order the baryon- 
to-photon ratio r. Numerically integrating over the whole /i-era, we estimate the effect of 
adiabatic cooling as 

f°° dz ( z \ 5 ^ 2 
fi = - -. -e~W 2.1 x r ~ -1.6 x 1(T 9 , (3.12) 

Jz f + 

where rather than a sharp cutoff at high redshift we use the exponential suppression at z\ 



obtained from the analytical solution of 26 ] . Notice that taking into account baryon temper- 
ature is essential in order to capture this effect. As we mentioned before, baryon temperature 
is a small effect and it is typically negligible unless all other effects cancel precisely. This is 
exactly what happens here: neglecting baryon temperature T oc a -1 , and so n 7 oc a~ 3 . Then 
the conservation of photon number is automatically satisfied by the solution even without 
enforcing N§ tft = and no /x-distortion is created. At first and second order waves propagate 
on top of the homogeneous background. Neglecting dissipation, each fluid element expands 
and contracts adiabatically, moving along the same adiabat as the background. Therefore a 
discussion analogous to the one above holds. Neglecting baryon temperature (or equivalently 
terms of order n&/n 7 ) no spectral distortion is created. In order to understand this from a 
slightly different point of view, let us discuss more in detail the conservation of entropy. 

3.2 Entropy considerations 

In this subsection we use again the argument presented in the introduction (|3.2|). Neglecting 



viscous corrections, entropy is covariantly conserved. We can therefore define 15 an entropy- 
density current = su^ 1 , where s = —u^s^ 1 is the rest-frame entropy density. This is 

15 A definition of the entropy density current that accounts for dissipative effect is in given in ( |C.13| ) 



- 21 - 



conserved during the adiabatic evolution and satisfies = 0. Both the entropy of baryons 
and photons is proportional to their number density 16 , hence up to terms of order n^/n^, 
discussed in the previous subsection, we can approximate s M ~ s 7 . Using Ts 7 = /? 7 +p 7 — //th?i 7 
one finds 

s 7 = \a R f z (1 + 39 r + 30 n(2) + 3Q 2 r + u° - A sfJ ) , (3.13) 



where A s = 135C(3)/(27r ). Combining this expression with (|3.5| ) one finds 

d *(— ) =£^{A n -A s )~l$(i. (3.14) 
Vn 7 / 3 6 R 

Using the conservation of entropy and photon number = = 0, we then find /i = 0, 
which result is valid up to correction of order fi^jn^. This result agrees with the discussion in 
the last subsection. As we will now see, once entropy is not conserved, things become more 
interesting. 

3.3 Viscous fluid 

Let us include the effect of dissipation during the //-era. The expressions ( |2.52| ) for the 



dissipative corrections obtained in subsection |2^2| are still valid, but we need to compute the 



heat conduction for the photon number conservation Xb- We leave the details of the derivation 
to appendix [D], and describes shortly the result. In the presence of dissipation, the velocity of 
energy transport differs from the one of number transport. Physically this can be understood 
as follows. Imagine two regions of the fluid with a different temperature. Suppose the same 
number of particles diffuse from the hotter to the colder region and viceversa. Then there is 
not net number transport, but there is a non- vanishing energy transport since particles coming 
from the hotter region carry in average more energy. There is another way of thinking about 
heat conduction. Consider a (very weakly interacting) gas of photons in a box. Consider 
the initial configurations in which the left hand side of the box has a certain temperature 
and number density different from those in the the right-hand side of the box. Imagine to 
choose temperatures and number densities such that energy density and pressure are exactly 
the same on the two sides of the box. For an ideal fluid, this configuration does not evolve. 
Conversely, once heat conduction is taken into account, a heat and number flow turns on and 
evolves the system towards equilibrium. 

The final result (see appendix for the heat conduction of photons is 

4 _ 2-7T 4 R 2 

^ r = 3^45037y' (3 - 15) 

where the apparent divergence as fl — > is fictitious since Ji cancels in the expression for 



^This can be checked using Ts x = p x + p x — ^ b n x . The baryon chemical potential is given in (2.58) 
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3.3.1 Zeroth order 

At zeroth order in perturbations, the only new effect with respect to the perfect fluid case of 



subsection \5A\ is given by the bulk viscosity, which is the only non- vanishing viscous coefficient 
at this order. Some distortion proportional to £ is generated, but as we discussed, this effect 
is very small, of order r 2 . In passing, this means that the background evolution is very close 
to adiabatic, i.e. no entropy is generated. In the following, we will neglect bulk viscosity. 
Since the background is still homogenous and isotropic, the relations Q2.57 ) are still valid. 



3.3.2 First order 



During the //-era, there is an extra equation in addition to those discussed in subsections |2.1.2 
and |2.2.2l JV^/j + AiV^ = 0. Using the results of appendix |] and and Q one finds 

405C(3) 2 
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AN* 



R 2 



7T 



6tt 4 C(3) 



-/' 



7 (l + i?) 2 
The resulting conservation equation is 



t 7 [d*Qr + &6n - u l H] 



-AA5fi + 4-^-u - 35, 



R 2 



7 



5fi + <L 



(3.16) 
(3.17) 
(3.18) 

(3.19) 



For modes well inside the horizon, as we did previously, we can neglect gravity perturbations, 
solve for u, u and 5b and obtain a system of a one first- and one second-order differential 
equations for Sfj, and <L. Notice thought that the longest mode inside the horizon during the 
/x-era is I ~ 1400 entering around Zf ~ 5 x 10 4 . We will discuss modes outside of the horizon 
shortly. By keeping only terms that are at most linear in either 5/j, or t~, inside the horizon 
one finds 

s, = - ,J n '3 3 L,,,f.\, (^-mL) (3.20) 
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3H5, 



(3.21) 



OR + 1) 2 

and for <5 7 the same equation as before the //-era, ( 2. 70]) , up to terms of order t^H and q/ (aH). 



Since the solution for <5 7 is again ( 2.71 ) at leading order, the above equation tells us how the 
oscillations of the fluid source short scale (I > 1400) perturbations in /i at linear order. 

Long perturbations in fi, i.e. I <C 1400, are of more direct observational interest. Since 
these modes were well outside the horizon during the /i-era, we can compute their evolution 
by neglecting derivatives. Since all viscous corrections (apart from the bulk viscosity which 
is extremely small) come with a higher number of derivatives, they are all negligible. Then 
the dynamics is the same adiabatic evolution that we studied in subsection |3.l| . Up to terms 
suppressed by the photon to baryon ratio, 5fi = 0. 
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3.3.3 Second order 



One way to derive an equation for the generation of /U-distortion at second order was presented 
in the introduction, subsection using the non-conservation of the entropy caused by 
viscous corrections. Here we present an alternative derivation that does not make use of any 
entropy consideration. 



As discussed at the end of subsection |2.2.3| , the leading sources for the second order 
perturbations arise when the time derivatives hit the dissipation factors in {uiU 1 ). Terms of 
order R are, on the one hand, small corrections during the ^-era and can be neglected. On 



the other hand, according to the discussion of subsection ^2, we know that up to viscous 
corrections, the system evolves adiabatically and (fi) = 0. The terms of order R in ( p. 74 ) 



are not viscous correction, e.g. since they survive as i 7 — > 0, and therefore can not source 
(/i). We verified this explicitly to leading and first subleading 17 order in aH/q, where q is the 
wave number of the short-scale dissipating perturbations. This is important since it ensures 
that all the energy stored in primordial perturbations is conserved during their evolution and 
is available at the time the perturbations reach the dissipation scale. Prom now on we set 
R = and neglect baryon conservation. 

Instead of taking the average as we did in subsection 2.2.3] and 2.1.3] , we work with the 
full energy and photon number conservation equations at second order 



(4ff + $2$ 

(3H + d t ) nL + cw. 



'(2) 



(3.22) 
(3.23) 



where the viscous corrections in these formulae can be omitted since they are higher order 
in derivatives. On the other hand, as we will see shortly, it is essential to account for the 
dissipation damping in the solution of the first order perturbations. To simplify the above 
equations we used a 2 HTl 
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-5n~,u l + u 
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we can solve ( p. 22 ) for ujL and substitute it into ( 3.23 ). The result is 

d t (4<5n 7 ( 2 ) - 3<5 7 ( 2 )) = (<5 7 u*) + 2d t (uiU 1 ) . 



(3.24) 
(3.25) 
(3.26) 

(3.27) 
(3.28) 



17 In order to have a cancellation at this order one should include in ( [2.39| ) a term subleading in aH/q, 
coming from (2.35) when the time derivative on the right hand side acts on the amplitude of <5 7 rather than 
on the cosine. 
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Using ( |3.5D and (3.6), we see that this is an equation just for fi: 



(i = (3A p - ±A n )~ l 2d t Uy + 3G 2 1} r ) + % {5. 



(3.29) 



As a check of this equation, let us take its expectation value, hence dropping the total deriva- 
tive: 

(A) = - 2 



AA n - 3A p 
0.93 dt Utfui + ZQ'^ 



d t ( (u l Ui + 39j 1)>r ) 



2 

(l),r> 



Using the relation between first order quantities @(i),r = ~ai we can rewrite (2.42) as 

(« , « i > + 3<e5 L)ir ) = |((^»p. 

The equation for \x simplifies to 



(3.30) 



(3.31) 
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This result agrees with the derivation in subsection 1.1 and can be interpreted as 

„2 



-1.40 x 



1 + w 



dm)) P , 



(3.32) 



(3.33) 



where we see that the origin of the factor of 3/4 found in |^Q| is the relativistic correction 
(1 + w)^ 1 to the energy of wave, which we discuss in appendix ||. For modes of observational 
interest one can drop the total derivative term in ( 3.29| ). Then integrating over time we find 



Ho(x, t f ) 



3A P - AA n L 



Uiu 1 + 3Q 2 (1 



(3.34) 



where the index "0" is a reminder that this is a monopole (the same in every direction) and [•];■ 
indicates the difference of its argument between the beginning of the /i-era at Zj ~ 2 x 10 6 and 
the end zj~5x 10 4 . Using ( |3.34l| ) one can compute the momentum space power spectrum at 
the end of the //-era (details are given in appendix [H|, see ( |H.4p ) in the limit of small momenta 
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(3.35) 
(3.36) 
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4. After the /i-era 

As the universe expands, after zt = 5 x 10 4 , kinetic equilibrium is lost since the momentum 
exchanged between photons and electrons is of order T/m e <C 1. From this moment onword, 
thermodynamic equilibrium is lost and no additional /i-distortion can be generated. In this 
section, we provide a transfer function to relate the /i-distortion at the end of the //-era to the 
one we can measure at late times 18 in the CMB. It turns out that the Boltzmann equation is 
particularly simple, so we abandon the hydrodynamic equations for this section and switch 
to the kinetic description. 

We consider the Boltzmann equation at linear order in perturbations, neglecting double 
Compton scattering and Bremssstrahlung in the collision term. We make a Bose-Einstein 
ansatz for the photon occupation number, with constant fj,, which is known to be a good 
description for not so low frequencies v/T 3> 10~ 2 . We neglect polarization and multiples 
above the dipole, since these are suppressed in the tight coupling regime. In general one has 
/j,(t,x,p,p) and our conventions for the decomposition in spherical moments are reviewed in 
appendix |F[ Going to Fourier space and introducing the cosine between the wavenumber q 
and the photon momentum (or direction of observation) p, x = P ■ k, we find 

m' + nxv = J-(IM}-H), 

where a prime denotes derivatives with respect to conformal time dr = adt. Notice that 
this is very similar to the analogous equation for temperature perturbations, with the major 
difference being that there is no velocity sourcing a dipole for /i. Following this analogy, we 
solve (f4.1[) in two steps. First we use the tight coupling expansion to find a solution in the 
epoch between the end of the /i-era and the last scattering surface, when the photon mean 
free path i 7 is the shortest distance in the problem. This gives us some /i monopole /j,q and 
dipole fjL\ at the last scattering surface. Second, using the integral along the line of sight we 
solve fl4.ip in the epoch between last scattering and arrival on the earth in the present day. 
This free streaming evolution relates fiQ and \x\ to the observed multipole at late times. 

4.1 Tight coupling 

Projecting on the zeroth and first Legendre polynomials, we derive the coupled equations for 
the \i monopole /j,q and dipole \i\ 

^0 + ^1 = 0, (4.2) 
fJ-i ~ 3M0 = -— W • (4.3) 

As usual, by projecting on higher Legendre polynomials, one can verify that in the tight 
coupling limit i 7 —> 0, the higher multiples are suppressed by powers of qt^ and hence can 

18 In jl3| this transfer function was taken to be just the projection of the monopole at the end of the /i-era 
onto higher multipoles in the late time CMB sky. As we show here, this is a very good approximation for large 
scales, I < 1000. 



- 26 - 



be neglected. We algebraically solve the above equations for fj,\ and and plug the solution 
back into the time derivative of the first equation. We find 



Mo +M07- + T^o = 0. 



(4.4) 



The two WKB solutions can be expanded in qt^ <C 1. The one that decays more slowly is 



fj, (q,t) = /x (<?,*/)exp 
= Mo(g,*/)exp 
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try dij 
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(4.5) 



tf 



where no(q,tf) is the value of the fiQ monopole at the end of the //-era, ( 3.34 ), and we intro- 
duced the decay rate and the /i-dissipation scale q^D, both functions of time. It is useful 
to estimate q^jj since it is the smallest scale on which we can observe some inhomogeneity in 
fi. We can use 
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2 J ' rrih 
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(4.6) 
(4.7) 



where the Helium fraction Y = 0.23 enters because of the relation between electron and 
baryon number density n e = (1 + Y/2) n^, ax — 66 fm 2 is the Thomson cross section, Hq = 
2.4 x 10 -4 Mpc -1 is the Hubble constant nowadays from it follows that today's critical density 
is p crit ~ (2.5 x l(T 3 eV) 4 . Then 19 



00 dz{l + z) 
1100 7 



H 

0.084 x Mpc -1 



-1/2 



(4.8) 
(4.9) 



corresponding to a multiple I^d — 1200, which is comparable with the damping scale of 
temperature anisotropies. 



4.2 Free streaming 

As free electrons combine with protons to form neutral hydrogen around z ~ 1100, the mean 
free path of photons grows larger and larger and the tight coupling approximation breaks 
down. In fact, most of the photons we observe in the CMB last scattered very close to 
z ~ 1100, i.e. have been free steaming since recombination. It is well known how to evolve 

19 The integral is supported at late times, so the upper bound of the integral is irrelevant as long as it is 
much larger than the redshift at last scattering zlls — 1100 
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temperature anisotropies during free streaming, and very similar techniques can be used to 
solve (14. IT). First we rewrite it as 



d_ 



K T > 9> X) exp (iqxr + J df J-^j = ^"M 7 ") q) exp (-iqxr - J df . (4.10) 



The line-of-sight solution is 

Kr, q) = I dfg(f, r)/i (f, q)e^^ , (4.11) 



TO 

where tq is some very early time and we have defined the visibility function (the same as for 
temperature anisotropies) 

We can project on the various multipoles by multiplying this solution times i l Pi(x) an d 
integrating over \- Using the identity 



J ^fl(x)e' : « <f - T) = 31 [q (r - f )] , (4.13) 



we find 



V>l(T,q) = J dfg(f,T)no(f,q)ji[q(T-f)] , (4.14) 
= lMi(g,r f )Mi(q,T) , (4.15) 



where in the second line we used the result of the tight coupling analysis (|4.5|) and introduced 
the transfer function 



Mi(q,r) 



J dfg(f, t) exp (- J ^ at^dnA ji [q (r - f)] (4.16) 



This integral can in principle be computed numerically with the exact visibility function. 
Analogously to what happens for temperature perturbations, we can try to consider small /, 
say I -C 1000 and approximate g as a delta function around the last scattering surface tlss, 
where hq and ji vary slowly. Then 

Vi{T,q) ^ IM){rf,q)ji [q(r - t LS s)} (tentative). (4.17) 

This is the transfer function used in [|l3|] and it is an excellent approximation for large scales 
for the goal of computing [iT correlations. For fi^i correlations, things are more complicated. 
In our formalism the late time angular power spectrum is 

Ct" = J 7^ P Mo(T/,?)|MKg,r)| 2 . (4.18) 
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where P^ijf^q) is the momentum-space three-dimensional power spectrum of po at the end 
of the //-era. The transfer function M;(g, Ti ate ) is effectively zero for modes q S> q^D{ z LLs)-, 
with zlss — 1100 the redshift of the last scattering surface. Hence we are only interested 
in the q — > limit of P^ , which is computed in appendix || and given in (H.4). In this 
approximation 

Cr = PM J ^|M^,r)| 2 . (4.19) 
If we approximated the visibility function as a delta function, we would get 

Ct" = Pvoirt) I ^ji{qr L ) 2 e-^l^ss? . (420) 

This integral is supported on large values of q, and cut off by the exponential diffusion damp- 
ing. But the thickness of the last scattering surface suppresses the small scale power spectrum 
by an amount comparable with the diffusion damping and so it can never be neglected, not 
even on large scales. An analytical expression can nevertheless be found by employing the flat 
sky approximation and approximating the visibility function as a Gaussian of conformal-time 
width cr 2 . Then one finds 

QD{Tf) r 2 L 

where l^o ~ fLlixD and q^D~ 2 = q~p + a 2 .. A few comments on this result are in order. The 
small scale pertubations present in two patches of the sky, separated by a large angular scale, 
I <C 1000, are independent random variables, so their correlation is just white noise. In fact, 
the p power spectrum is I independent for small I, as expected for white noise. 
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A. Conversions between single fluid and two fluid descriptions 

The photon-baryon-electron plasma froms a single fluid with p = p y + and p = p 7 = p 7 /3. 
We define 20 

P 1 2 _ ( d P\ 1 



W -p 3 + 4.R' Cs ~{dp) s ^ Nb 3(1 + R) (AA) 



20 In principle the derivative defining the speed of sound should be performed along an adiabatic trans- 
formation, but since all the processes we consider are reversible, this is the same as constant total entropy 
S. 
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where S is the total entropy and Nb the total number of baryons. The fastest way to compute 
(? s is to notice that the background evolves adiabatically conserving the baryon number. 
Therefore one can use the chain rule to write (? s = p/p. This gives the correct result. To 
make more transparent the role played by the conservation of entropy and baryon number, 
in the following we provide an alternative (longer) derivation. 

There are six relevant thermodynamical variables: S, N b , p, V, T and U. By the defini- 
tion of (? s we are instructed to keep the first two fixed and take a derivative in the direction 
p = U /V . To compute this derivative we need other three independent equations that relate 
the extra three variables. The first equation is given by the integrated form of the second law 
of thermodynamics 

TS = U + P V - p b N b . (A.2) 

A second an third equations are given by explicit expression for the pressure p and entropy 
density s = S/U. Neglecting baryon temperature, one has 

S = Il5 T3 ' P=l(p-^n b ), (A.3) 

where we have used the energy density p = U/V and baryon number density n b = N b /V. The 
equation for s comes just from the black body formula, while the one for p uses p = p 7 + p b 
and from (2.58) p b ~ m b n b ~ p b n b . Notice that this relation implies that the photon-electron- 



baryon plasma is not a barotropic fluid, i.e. p 7^ p(p)- Using these three equations one finds 

{d p log s) SjAr = (d p log n b ) S N = - {d p log V) SjN , (A.4) 

which can be solved for (dT / dp) <j N . Solving (|A]| and flAj) for n b as function of p and T 
and computing the appropriate derivatives, one derives the well-known result 

where as usual we have defined the baryon loading R = 3pfe/(4p 7 ) for which R = HR. 
Some useful relations are 

V-uzh. §4 = -T#7?' < A ' 7 > 

ur Hcj 1 + R 



B. Initial conditions 

In this appendix we derive the initial conditions for the amplitude of the pressure waves in 
the fluid. We are interested in radiation domination so we neglect dark matter. We also 
neglect neutrinos which will be discussed elsewhere. Their effect is about 10%. So we are left 
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with a single photon-baryon-electron fluid (in the tight coupling regime) coupled to gravity. 
Following plj we perform the computation using a P(X) Lagrangian which is equivalent to 
a single perfect fluid. We work in comoving gauge, i.e. u = for u the velocity potential and 
neglect tensor and vector perturbations. Notice that in this gauge 5p 7^ 0. After solving the 
constraints from gravity one gets the simple second order action 

„2 



S2 



dtd A xa 6 ^ 



% {dizf 



where 1Z stands for curvature perturbations on comoving slices 21 , 
tion is 

A 

Tl=- + HSu, 



(B.l) 

Its gauge invariant defini- 
(B.2) 



where as defined 22 in ( 2.22| ), 5gij = a 2 (A5{j + didjB). During radiation domination e = 2 



and for the photon-baryon-electron plasma c~ 2 = 3(1 + R), with the baryon loading R = 
3p6/(4p 7 ). Much before matter-radiation equality R >C 1 so we neglect it in the following. 
The equation of motion for 7Z in Fourier space is 

U" + 2aHK' + q 2 c 2 s K = , (B.3) 

with the conformal time r = J dt/a = l/(aH) = 2t/a. Notice that r > and runs from 
zero to positive infinity. The solution of ( B.3 ) asymptoting a constant in the far past 
(r, a — > 0) is 



n 



qc s 

a 



1 



^ (0) sin 

qc s 2t 



dt 



(B.4) 
(B.5) 



Notice that with the definition ( |B.2| ) , this result is valid in any gauge, both inside and outside 
of the horizon, during radiation domination. If we consider modes well inside the horizon 
q/(aH) 3> 1, synchronous and Newtonian gauges agree. In either of these two gauges, the 
gravitational potential is very small on scales smaller than the horizon so we have 1Z ~ H5u. 
Again neglecting gravity perturbations, the conservation equation for the fluid is 

„2 



5 = (1 + w)—^5u . 



Integrating we find 



1 + w 



TZ° cos 



-c, dt 



4ft° cos 



-Co dt 



(B.6) 

(B-7) 

(B.8) 
(B.9) 



which agrees with (6.4.11) of p| and (121) of p| 



21 This quantity is called £ in Q 
22 This is the notation of flj^. 



- 31 - 



C. Derivation of the viscous parameters 



For some conserved charges Na with chemical potentials p,^, the second law (sum over 
repeated index A = 0, . . . , q) 

TdS = dU + pdV - p? h dN A , (C.l) 
can be rewritten in terms of densities s = S/V, p = U /V and ha = Na/V. Then 

Tds = dp- pth dn A + ^-{p + P- pth n A ~ Ts) . (C.2) 

The volume is so far unspecified. It is convenient to define the volume such that it contains 
a fixed number of say Nq particles. Then dNo = and hence VdnQ = —riQdV. We introduce 
the entropy per particle of type "0", a = s/uq, drop the index "0" so that its number density 
will be simply called n, and define a = 1, . . . , q. Then one can rewrite the second law as 



Tdc = d(P-)+ P ( l -)-pl h d( T ± 
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dp 
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nd^a 
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p + p 
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p + p 



n 



dn 



dpji - pt h n d. 



n,, 



n 



(C.3) 
(C.4) 

(C.5) 



Leading viscous coefficients 

It is convenient to we work in flat space and covariantize the final result to get the general 
relativistic expression. Let us consider the conservation laws 

= o = d^ A , (c.6) 

for every A = 0, . . . , q and define as usual 



= (p + p) u»u u + if v p + AT^ V 
iV£ = n A u" + AN* , 



(C.7) 
(C.8) 



where u^u,, 



1 and we work in flat space. As discussed in |2.1| , p = —u fl u u T fJ,u and 

|15[ . By multiplying the conservation 
of T^ u by u u and using the conservation of Nq (remember we will be dropping the label "0"), 
one obtains 



n = —u^N^- and we can choose u l such that u^AT^' 



p + p 



dfj,n 



-AT^d„u u + ?±^d^AN» 



where we used u^d u AT^ v = —AT^dyU^. Using this in (I 
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one finds 



nw 



l d^a = -T-'AT^d^ + P±^d^AN» , 
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u»d„— . 
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(C.9) 



(CIO) 
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We would like to rewrite this expression as 



<9 M s M = positive terms , (Cll) 

where s M is the entropy density current (which does not make reference to any conserved 
charge). In this way, the entropy can only grow for an arbitrary flow. The difficulty is finding 
the right velocity to define the current s^. We adopt the following strategy. Write some 
ansatz for and solve for the free parameters in the ansatz such that on the right hand side 



of ( Cll ) depends on AT and AN a but not on their derivatives. Once this is done, we ask 
that these terms form perfect squares. This determine AT and AN a- 
We start with the ansatz 

S M = sn M + f AAN » ^ ( C . 12 ) 

for some arbitrary functions /a and the sum over A is implicit. We could not have added any 
other term, since u (i AT' 11 ' vanishes because of our choice of u^, and using derivatives gives 
subleading terms in the hydrodynamic expansion. Physically this ansatz say that the entropy 
density flow with the velocity of energy transport plus some corrections accounting for the 
difference between and the velocity of particle transport. Plugging ( |C.12 ) into ( C.1C ) one 
finds 



V = -T^AT^d^ + AN^fA (C.13) 

where we remind the reader that n = ua=o- Choosing Ja = — /i^/T gets rid of the second 
line and with it of all derivates acting on AT and AN. So the final answer is 

s m = su m_^AA^, (C.14) 
= -T^AT^d^ - AN»d^. (C.15) 



The first term on the right hand side of (C.15) captures bulk and shear viscosity corrections 



to T^ v . For these terms the discussion is the same as in |{L7| , so we do not repeat it here. 
Let us focus instead on AiV. A sufficient condition for the corrections to be always positive 
is AN^ oc d^fx^/T) . From now on the discussion is the same for every A, so we drop the 
specie index. From the definition of n, one must have u^AN^ 1 = 0. To enforce this we can use 
the projector H^ v = 7]^ v + u^u v in AA^ cx H^du^fj^/T). The constant of proportionality 
is a matter of convention and dimensional analysis. We use the convention of ( 2.52Q . 



D. Derivation of the heat conduction 

In this appendix we derive the formula for the heat conduction in the conservation of the 
number of photons x-y an d baryons Xb- The latter has been known for a long time [], while 
the former, to the best of our knowledge, has not appeared yet in the literature. 
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Baryon heat conduction 

Let us start with the Boltzmann equations for photons and baryons. Notice that we can carry 
on this computation in fiat space and use the result in an expanding background with small 
perturbations. Since in the simplest models the spectral distortion is expected to be much 
smaller than the other perturbations, we will neglect it in this computation. The first two 
momenta of the two relevant Boltzmann equations are [25| 

e + g9i = 0, (D.l) 

Vb = ^ K + 3i0i] , (D.3) 
5n b + iqv b = , (D.4) 

where, using the notation of 25] r' = —t^ 1 and we have defined the baryon velocity potential 
by v\ = v b q l /\q\. The idea is now to find the difference Au b between the velocity appearing 
in the conservation of baryon number and the one appearing in the conservation of the total 
energy density u. So we write 

5 = q 2 (l + w)u, 5n b = q 2 (u + Au b ) , (D.5) 



where w is given in (A.l), 



(D.6) 



denotes perturbations to the total energy density and 5n b is dimensionless. Using ( p.4j ) and 
( p.5p one finds 

q 2 Au b = - (iqv b + q 2 u) (D.7) 
One can solve ( p.3| ) for v b in the tight coupling expansion. At first order the result is 

iq Vb = SqQ, + + O (f~ 2 ) . (D.8) 

Using the other Boltzmann equations and the definition of u one can similarly find 

q 2 u = -3q& 1 - 0§±r + O {f~ 2 ) , (D.9) 

■ ge q R 2 e 

0l " WThT) ~ 3(1 + R) 2 r ■ (D - 10) 



Hence finally 
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One can recast this using the notation of pq] , as in ( 2,52 ) , provided that 



XbT = ^p 7 *7 , (D.12) 

which agrees with [jl7|]. 
Photon heat conduction 

In order to compute the heat conduction x-y appearing in the conservation of photon number, 
we could use the same strategy as above. On the other hand, there is a shortcut that leads 
to the same result. Let us start by writing the transport of energy density in two different 
but equivalent ways 

(p + p) u = (p 7 + p 7 ) (u + Au 7 ) + {pb + Pb) (u + Aiifc) . (D.13) 

By using the results of the previous section we can then find 

R 2 

A-u 7 = -RAu b = (1 + jR)2f Qo • (D.14) 

During radiation domination, when R <C 1, the velocity of energy transport is closer to the 
velocity of photon rather than baryon number transport. This is intuitive since the photons 
carry most of the energy. 

Again one can recast this using the notation of fllql , as in ( 2.52j ). The resulting heat con- 



duction appearing in the conservation of photon number, caused by the presence of baryons, 
is 

4 2tt 4 R 2 

4 R 2 
~ -p 7 i 7 3.6 x — (D.16) 

where the apparent divergence as ft — > is fictitious since p cancels in the expression for AN 11 
or A«~. 



E. The energy of a wave from fluid dynamics 

In this appendix we derive a relativistic formula for the energy density of a pressure wave in 
a perfect fluid. For simplicity we work in flat, unperturbed space. To begin we have to give 
a precise definition of the energy of wave up to second order in perturbations. For a perfect 
fluid we define 

E w = (AT 00 )p , 

AT 00 = (T 00 [(p) = p av y]) - (T 00 [(p) = p av y = 0]) (E.l) 
= (T 00 [(p) = Pm J])- P av, (E.2) 
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for some average energy density p av . Here we have introduced the average over a wave period 
In words, the energy of a wave is the energy of the fluid with the wave (v? ^ 0) minus 
the energy of the fluid without the wave (u l = 0) with the same average rest-frame energy 
density p av . Expanding in small perturbations around a homogeneous and isotropic solution 
(potentially time dependent) we define 



P = P 



l + 5 + 5 {2) ) = p + Sp + 5p (2 ), 



+ m + u 



(2) 



so that at second order 



T 00 = p (l + 5 + # 2 >) + (p + p)u l u t 



Using the definition (E.l) one finds 



En 



{AT 



00 \ 



{{p + p)u % Ui) p . 



(E.3) 
(E.4) 

(E.5) 
(E.6) 



In order to rewrite this expression in terms of density perturbations, let us consider the 
dynamics. There are four equations T£ u = for five variables p, Ui and p. An equation of 
state will be needed to close the system. For a perfect fluid fl2.4| ), energy and momentum 
conservation can be written as 



dt [(p + p)(l + 



0, 



dt 



(p + p)uiy/l + |li| 2 + 8j [(p + p)u j Ui + Sijp] = , 



(E.7) 
(E.8) 



where \u\ = u l Ui. We solve these non-linear equations in perturbation theory around a homo- 
geneous and constant background. For this purpose it is convenient to take the time derivative 



of |E.7| and solve for dtdi 



pui^Jl + 



u 



. Then one substitute this into the divergence of E.5. 



At linear order and going to Fourier space, the result is 

dp + q 2 5p = . 

From ( |E.S| ) focusing on scalar degrees of freedom one finds 

(p + p)q 2 u = 5p. 

Barotropic fluid 



(E.9) 



(E.10) 



If we assume that the fluid is barotropic, then the equation of state takes the form p = p(p). 
This in particular means that 



Sp = = c 2 s 5p, 



(E.ll) 
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where we used the definition of the speed of sound (Ad). Notice that if the fluid had not 
been barotripic there would have been additional terms not dependent on c 2 . The photon- 
baryon-electron plasma is not in general barotropic, since (neglecting baryon temperature) 
p = p 7 (4/3)(l + R) and p = p-y = p-y/3. On the other hand, it can be approximated as 
a barotropic fluid in at least two case. At very early times, neglecting terms of order the 
baryon loading R. During an adiabatic evolution, since then the system moves along a one 
dimensional adiabat, which can be parameterized by p. Hence the following formulae will be 
a good description of the photon-baryon-electron plasma in these limits. 



Substituting ( E.ll ) into ( [E.9D one finds oscillatoric solutions, which in the WKB approx- 
imation take the form 

5p = A cos (^j qc s dt' + phase^ . (E.12) 

Hence, always at leading order in the WKB expansion (i.e. we neglect the possible time 
dependence of the amplitude ^4) 

(6p 2 ) p = q 2 c 2 s (6p 2 ) p . (E.13) 



Using ( |E.10D and (Je7T|) we find 



c 2 

E w = {{p + p)u l Ui ) p = ^f^(Sp 2 ) p (E.14) 

P + P 



1 + w 



(E.15) 



where in the last line we used w = p/p. This result agrees with the non-relativistic formula 
for p <C p. For radiation c 2 = w = 1/3, and hence 

E w = (3/A)pc 2 (5 2 ) = \p(S 2 ). (E.16) 
F. The energy of a wave from kinetic theory 

According to p5|], unlike other folks physicists have always one thing in mind: the Boltzmann 
equation. Let us therefore verify the hydrodynamic result ( |E.16| ) using kinetic theory. We 
start introducing our conventions and move on in the next subsection to discuss the formula 
for the energy of the wave. The energy-momentum tensor in an FLRW universe for a gas of 
bosons is 

TV (3 t) = - ! l^tl_ \ rp i) 

K,) a 3 J (2tt) 3 p° eP °/kT(x,t,p) _!> ^- L > 

where p°p° = piPia~ 2 = p 2 a~ 2 and pi = Pi/p. It is common to denote small temperature 
inhomogeneities by 

T(x,t,p) = T[l + e(x,t,p)] , (F.2) 
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where T is a dimensionful constant and for the moment we do not allow for distortion, i.e. it 
does not dependend on p, but just on p. From now on we omit to write the time dependence. 
We can decompose the Fourier transform 



Q(k,p) = J , 



d 6 xe ikiX ' Q(x,p) , 



in multipoles 



e(k,p) = G(k, x ) = E ( 2/ + x ) (-*)' p i(x)&i(k) , 

1=0 

&i(k) = i l J 1 ^e(k,x)Pi(x) = * 1 J ^@(k,p)Pi(oos(9 p )), 
where x = kipi/(kp) = cos (9k P ) and the Legendre polynomials satisfy 



In the following we will also use 



d 4 k 

da 



47T 



^e(x,p)p ; (cos(%)). 



With these definitions and conventions one finds 

72 ( 



r oo (f) = pj d -^[l + 4Q(x,p) + 66 2 (x, p)] , 



(F.3) 



(F.4) 
(F.5) 



(F.6) 



(F.7) 
(F.8) 



(F.9) 



where p = fg- (fcT r ) . Decomposing in multipoles and using the statistical homogeneity of the 
perturbations one obtains 



(T 00 ) = P 



i + 4(e (x)) + 6 J — £(2/ + i)<|e^)| 2 >' 



(F.10) 



where 



(2^) 3 ,5 3 (k+k>) (le.wi 2 )' = mmik')} 



(F.ll) 



and we used the reality of Oq^O?) to substitute 0o,r(— = ©o ir (&). 
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The energy of a wave 

In the tight coupling regime at every point we can go in a boosted reference frame which we 
denote with the label "r" (since this is the frame in which the baryons are at rest) in which 

T(x,t,p) = T r [1 + e , r (x)] , (F.12) 

i.e. there is only a monopole. Since we assume statistical homogeneity (@Q tr (x)) does not 
depend on x and can hence be absorbed into T r . So without lost of generalities we impose 
(@o t r(x)) = 0. T 00 (x) in this comoving frame is p(x). For small inhomogeneities we have 

T r 00 (x) = p[x) = p[l + 49 , r (x) + 6Ql r (x) + O (0 3 )] , 



and therefore 



( P (x)) = p[i + 6(el r (x)) + o{e 3 )] 



1 + 6 



d?k 
(2tt)3 



(\e , r (k)\ 2 y + o(Q 



(F.13) 

(F.14) 
(F.15) 



Notice that p is defined as the energy density measured by an observer that sees the fluid 
locally at rest. Let us now review how temperature changes as we go to a boosted reference 
frame. If we make a boost with velocity v = Vi, say from reference frame S to S' , the 
directional temperature changes from T(p) to T{p') according to (see e.g. f27j) 



T{P) = T{P) VT^ 



with 



v ■ p 



v ■ p + v 



(F.16) 



(F.17) 



1 + v ■ p' 

So going from the rest frame ( |F.12| ) to a boosted frame, one finds up to quadratic order in 
velocity 



Q(x,p) = Qo,r(x) + v- p+(v-p) z 



+ 0{v A 



(F.18) 



This tells us that the boost creates both a dipole and a quadrupole and modifies the monopole. 
More specifically, let us focus on an irrotational velocity field Vi(x) = diu(x), so that Vi(k) = 
iku(k). We find 

d 3 k 



e (aO 
ei(fc) 



Ait 



-Q(x, p) 



1 



(2vr)= 



3 ik ^ e (k) = e , r (x)--v 2 (x), 



6 



1 



-k u(k) 



If we substitute this result into flF.10 ), using (0o, r (x)) = 0, we find 



(T 00 ) = p\l + 



d 3 k 
(2tt) 3 



6(|e (A:)| 2 > / +(^-^)^(KA ; )| 2 )' 



(F.19) 
(F.20) 

(F-21) 
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Using the definition of appendix |E| and ( F.21 ) and ( F.14 ) we find 

d 3 k 4 



E„ 



(2vr) 3 3 



k\\u{k)\ 



2\i 



p-(ui(x)u l (x)} , 



(F.22) 
(F.23) 



which, for radiation, agrees with ( |E.6| ). 

G. Evolution of distortion in the cosmological frame 

In this appendix we re-derive an equation for the time evolution of /U in the cosmological 
reference frame, in which the fluid as a whole does not move. A parallel derivation using 



rest-frame quantities was given in subsection 1,1 . The ensemble-averaged conservation of the 
energy- momentum tensor at second order is given in ( 2.43j ) 



(dtT$ + 3#Tgj + a 2 HT { 



oo 

(2) 



-(2)/ 



0. 



(G.l) 



With some insight we have neglected the viscous term AT^ , which is subleading well inside 
the horizon. Approximating the plasma as a gas of photons, we can use T VL = a _2 T 00 . The 
conservation of photon number gives 



(d t N° + 3H N°) = . 
Let us expand in perturbation according to 



-.00 



fW(l + 5T™+5T™ 



N u = N u (1 + 5N^ + 5Nfo 
Then at second order we find the simple ensemble-averaged equations 

d t (6N° 2) )=d t (6T$)=0. 
Let us now use the expressions from kinetic theory 



T 00 (x) = (kT^ 



^ + 40(1)+ 4e (2) , o -^ + 6^ 



N°(x) = (kT ) : 



15 



;2C(3) 



71 



7T 1 



7T 



1 + 3G (1) + 36 (2)i0 - _ 



M + 3 



4vr ^) 



where the temperature in the cosmological frame is 

T(x,p) = T [l + 9(i) (x,p) + 8(2) (x,p) 
and we have used the notation 







(2),0W 



4 7T 



0(2) (X,P) 



(G.2) 

(G.3) 
(G.4) 

(G.5) 

(G.6) 
(G.7) 

(G.8) 
(G.9) 
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Substituting (|G4) and (p~7|) into (|G~5|) , one finds 

9 'W = - 2P^(3FI x29 '/§< 9 a>) < G10 > 

^3 [(l©(i),o| 2 )+3(|e (1)il | 2 )] (G.ll) 
= -1.40xi^), (G.12) 

where in the last line we used the wave solution to relate first order monopole and dipole 
9(i),i v3 = ©(i),o = (see e.g. This result agrees with the estimate 



dtfr) = -1.40^ ~ -1.40-^-^(4), 
for c 2 = w = 1/3. The time evolution of the temperature is 



7^3 [(l©(D,o| 2 )+3(|e (1)il | 2 )] (G.15) 



-1.14 x -9 t (^) . (G.16) 



The above derivation, performed in the cosmological frame, agrees with the one presented in 
subsection |1.1| performed in the rest frame. First, \i does not depend on the reference frame, 
and in fact ( |G,13 ) is identical to ( 3.32 ). Second, given that the temperature transforms as in 



( |F.16 ), we expect 



<e (2) ,o - e (2) , r ) = ~(v 2 ) = ~x \d t (5 2 ) . (G.i7) 

From 

(405C(3) - 2vr 6 ) d^m) + 9 (270C(3) 2 - tt 6 ) a t (6 2 1) r ) 



^(@(2), 



6 [vr 6 - 405C(3) 2 
and ( G.14| ) one can check that this is indeed the case. 



0.59d t (u l Ul ) - 2.3^ (GL _). (G.18) 



H. The /i power spectrum 

In this appendix we give a derivation of the three-dimensional power spectrum of [i at the 
end of the /i-era. This can be used together with the transfer function, which we discuss in 



subsection 4.2, to obtain the late time /-space power spectrum. As explained in subsection 



only the power spectrum (/j,(q) fj,(q)) in the limit q — > affects observations on angular 
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scales relevant for observation. Hence in the solution ( p.29|) , we can drop the last term since 
it is a total derivative. Then we can easily integrate over time and find 



fM)(x, tf) 



3A p - AA n 



Uiu 1 + 3G 



(!),»■ 



/ 



(H.l) 



where tt = t(zf) and [•];• indicates the difference of its argument between the beginning of the 
//-era at ~ 2 x 10 6 and the end Zf ~ 5 x 10 4 . We now use 6(i), r = &y/4, 4qu/a = 3<5 7 and 
the explicit expression for <5 7 in ( 2.71 ). Neglecting terms of order R, the Fourier transform of 
fj, is then 



p 

sin 
+ cos 



3A p -4A n J (2tt) 3 
*t k 



d 3 k 

9c 2 s TZ (k)TZ°(q- k) (1 - 0.268i^) 2 



a 
'/ k 



c.5 dtj sin ^ 
Ccdt' ) cos 



dt' e-( k2 ^~ k nin + 



(H.2) 
(H.3) 



a 



k\ 



-c.c e 



(k 2 +\g-k\ 2 ) q - 2 



We can use this expression to compute the momentum space two-point correlation function 
in the limit q — > 



2(1 -CJ.268^) 2 f d 3 k2iT 2 A^(k) 
3{3A p -4A n ) J 17 fee 



-2k 2 /ql 



For a primordial power spectrum close to scale invariance, the integral is manifestly supported 
on the largest values of k. For an exactly scale invariant primordial power spectrum of 
amplitude A^(k p ), this reduces to 



P po (q^0, Tf ) 



2(l- 0.268^) 2 f d 3 k 2vr 2 A 4 e (fc) 
3(3A p -4A n ) J 17 



A: 6 



~ 2 x 10 



-14 



2(1 - 0.268.R,. 



2vr 2 A 4 ,(fe p ) 



3 (3A p - 4A n ) 
4 2(1 -Q.268^) 2 



2n 2 AUk p ), 



q D {z f f 3{3A p -4A r 



(H.4) 
(H.5) 
(H.6) 



which we simply call P^ {rf). 
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